Transverse instabilities of multiple vortex chains 
in superconductor-ferromagnet bilayers 
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Using scanning tunneling microscopy and Ginzburg-Landau simulations, we explore vortex con- 
figurations in magnetically coupled NbSe2- Permalloy superconductor-ferromagnet bilayer. The 
Permalloy film with stripe domain structure induces periodic local magnetic induction in the su- 
perconductor, creating a series of pinning-antipinning channels for externally added magnetic flux 
quanta. Such laterally confined Abrikosov vortices form quasi- ID arrays (chains). The transitions 
between multichain states occur through propagation of kinks at the intermediate fields. At high 
fields we show that the system becomes non-linear due to a change in both the number of vortices 
and the confining potential. The longitudinal instabilities of the resulting vortex structures lead to 
vortices 'levitating' in the anti-pinning channels. 

PACS numbers: 74.78.Na, 74.25.Qt, 64.60.Cn 



Recently, much attention has been devoted to the stud- 
ies of hybrid systems comprising of two or more elements 
with complementary physical properties, with the moti- 
vation that the resulting hybrid has a superior perfor- 
mance as compared to its constituents. Superconduc- 
tor - ferromagnet hybrids exploit the interaction between 
the two 'reservoirs' of strongly correlated electrons result- 
ing in a wealth of new physical phenomena (for review 
see The fundamental premise has been to isolate 

and tailor the dominant interaction between two cor- 
related systems, as many overlapping interactions may 
conceal the discovery of new physical phenomena. Re- 
cent explorations focused solely on magnetic interaction 
between superconductor and ferromagnet have led to 
the discovery of numerous intriguing phenomena asso- 
ciated with magnetic pinning of Abrikosov vortices as 
well as with mesoscopic confinement of superconductiv- 

ity aaiaasn. 

Confinement of superconductivity to quasi- ID trans- 
lates into confinement of the superconducting vor- 
tices. Abrikosov vortex structures constrained to one- 
dimensional superconducting condensate has been ac- 
tively studied 0, Ell El EI El, EI El in order to ex- 
amine novel static and dynamic vortex phases. Magnetic 
imaging of the vortex structures in these systems have 
been challenging, in particular in the cases where effec- 
tive penetration depth of the superconductor (i.e. mag- 
netic signature of the Abrikosov vortex) is on the order 
of the intervortex spacing. On the other hand, scanning 
tunneling microscopy (STM) can be successfully used to 
image the distribution of the local order parameter on 
the surface and map the Abrikosov vortex distribution 
in higher magnetic fields. 

In this paper, motivated by the aforementioned results, 
we study the vortex configurations in magnetically cou- 
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FIG. 1: (Color online) (a) Fabrication of magnetically cou- 
pled superconductor-ferromagnet structure: on a single crys- 
tal NbSe2 a thin layer of insulating Si02 is evaporated fol- 
lowed by a \\im thick Permalloy. After cleaving the NbSe2 
down to submicron thickness, fresh atomically flat surface is 
exposed to the STM probe, (b) Magnetic force microscope 
images of the stripe domain pattern in Py at Oe and 2000 
Oe applied field perpendicular to the surface of the film at 
room temperature. Scan area is 20x20 /mi 2 . 



pled superconductor (NbSe2) /ferromagnet (Permalloy) 
bilayer. We take advantage of well-ordered magnetic 
state of the Permalloy (Py) film to induce ID vortex 
confinement in the adjacent superconductor. We com- 
bine STM and Ginzburg-Landau simulations to explore 
the magnetic interaction between the magnetic domain 
state in a ferromagnet and Abrikosov vortices in an ex- 
treme type-II superconductor. NbSe2 was chosen for its 
small coherence length (compared to the ferromagnetic 
domain size), negligible intrinsic pinning, and atomically 
flat and inert surface. The choice of thick Py film with 
well-ordered stripe domain pattern allows us to obtain a 
ID periodic potential modulation of the superconducting 
condensate. Earlier we have shown that such magnetic 
potential emanating from the ferromagnet can be used 
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to achieve domain wall superconductivity [8]. Deeper 
in the superconducting state, strong interaction of the 
Abrikosov vortex lattice with the periodic magnetic do- 
main structure leads to commensurability effects [171 ]. 

The bilayer for present study was prepared using the 
method elaborated in Ref. [lfj. 20 nm thick Si02 in- 
sulating thin film was sputtered using RF magnetron 
on a freshly cleaved high quality NbSe2 single crystals 
(T C =7.2K, RRR=40). This layer was covered with a 
D = 1 jira thick Py (Fe2oNi 80 ) over layer using DC mag- 
netron sputtering from a single target (Fig. []Ji). As- 
prepared sample was attached to a conducting substrate 
with the ferromagnetic layer down. The surface of NbSe2 
single crystal was cleaved successively until the color of 
the surface started changing. This is a first indication 
that the thickness of the NbSe2 top layer is within the 
wavelength of the visible light. By further cleaving the 
sample the areas of different color could be obtained 
that signify different thickness of the superconductor [18[ . 
The final thickness used in the STM studies in this work 
was c.a. 0.4 jiva and this value was verified by cross- 
sectional cut of the sample using focused ion beam etch- 
ing and scanning electron microscopy metrology. 

The thick Py film forms narrow magnetic stripe do- 
main pattern with a period comparable to the film thick- 
ness 0, Hj (Fig. [Do). The period of the domain struc- 
ture W is proportional to the thickness D of the Py film. 
The magnetic stripe domain structure is insensitive to 
magnetic field applied perpendicular to the film's sur- 
face up to 400 Oe, but forms checkerboard domains in 
higher fields (Fig. [TJo). In our theoretical calculations, 
we used the dependence W « 0.0109V^D cm, following 
from Ref. [2l| and estimated magnetization of the film 
M = 1620 G, exchange constant A = 2.665 • 10 -6 erg/cm 
and uniaxial anisotropy constant K = 2.63 • 10 5 erg/cm 3 . 

The theoretical simulations of the described structure 
were performed within Ginzburg-Landau (GL) theory. 
We solve self-consistently a set of mean field differen- 
tial equations for the order parameter ip and the vector 
potential A: 



(-iV- A) 2 ^ = 
-ft 2 V x V x A 



(1 



(i) 

(2) 



The latter is the Maxwell- Ampere equation with a cur- 
rent density j = — \^\ 2 A. GL constant k equals 
the ratio between the magnetic field penetration depth A 
and coherence length £, and is in general very large for 
NbSe2 samples. In above expressions, all distances are 
expressed in units of £(T = 0), the vector potential in 
</>o/27r£(0), and the order parameter in ^/—a/f3 with a, (3 
being the GL coefficients. For more details on the numer- 
ics and implementation of periodic boundary conditions 
we refer to Ref. [22[ . Note that virial theorem requires in- 
teger number of flux quanta in the simulation region with 
periodic boundary, and we accommodate that by fixing 



the width of simulation region to the period of magnetic 
domain structure W, and the length is fitted (at « 2W) 
according to the applied magnetic field. 

In STM measurements, the vortex images were ob- 
tained by spatially mapping the tunneling conductance 
value at the edge of superconducting energy gap at 4.2 K 
fl5| | . In Fig. [2j we show the local density of states (LDOS) 
map of 0.7 x 0.7 jim 2 area on NbSe2 surface. Fig. [2] (a) 
shows a slightly distorted vortex chain at -100 Oe. The 
effective magnetic field above each domain is a sum of 
the contributions of the applied magnetic field and the 
field due to the local magnetization of the Py. In the case 
when these two contributions are of the same sign the ef- 
fective fields are enhanced, and when they are of opposite 
sign the two contributions cancel each other. From our 
map it is obvious that the central part of the image is 
above a domain that has a "negative" polarity, i.e. the 
domain is collinear with the negative applied magnetic 
field. One expects that at 100 Oe the average vortex dis- 



tance should be a v 



0.489 fim, and 



instead, at -100 Oe the vortices form a quasi- ID chain 
state with average distance of ~ 0.210 /im. No vortices 
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FIG. 2: (Color online) (a) STM images of vortex configu- 
rations in the NbSe2 at 4.2K. Applied magnetic field values 
(in Oe) perpendicular to the surface of the superconductor 
are shown in the upper right corner. The white dotted lines 
show the underlying magnetic stripe domain boundaries. The 
scanning area is 0.7 x 0.7 /im 2 . (b) Cooper-pair density 
plots, calculated at -105 and -165 Oe (in comparison with 
(a), see marked area), obtained for D = 1 /am (plot size 
1.087 x 2.176 /im 2 and 1.087 x 2.193 /xm 2 , respectively). 
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FIG. 3: (Color online) (a) Same as Fig. [2] but for H = ±300 
Oe and ±400 Oe. (b) Cooper-pair density plots, calculated at 
-325 and -385 Oe, for D = 1 /xm (plot size 1.087 x 2.176 /mi 2 
and 1.087 x 2.193 /mi 2 , respectively). 



are observed at +100 Oe in the same area, being pinned 
at the adjacent domain. 

Further increase of the negative applied field leads to 
increase in vortex density and formation of chain states 
similar to ones predicted in Ref. Q and observed earlier 
in patterned superconducting NbSe2 [15]. The chains 
consist of periodically spaced vortices in the middle of 
the stripe domain. As magnetic field increases so does 
the vortex density, leading to smaller inter- vortex spacing 
within a single chain. This process is continuous until a 
critical point leading to geometrical ordering transition 
and formation of double chains (Fig. [2(a) at -200 Oe) and 
triple chains (Fig. [5(a) at -300 Oe). The similarities in 
vortex ordering and geometrical phase transitions with 
the case of strong ID confinement potential 0, [l5] end 
at this point. Further, with the help of simulations, we 
focus on the subtle but important effects that are specific 
to the nature of our confining potential. 

We take advantage of the numerical simulations to re- 
construct the vortex phase diagram of the magnetically 
coupled S/F bilayer system. Numerical solution of Eqs. 
([T][2]) minimizes the Gibbs free energy of the sample, and 
Fig. [H shows the main outline of the ground-state phase 
diagram as a function of thickness of the Py film. Since 
this is a 2D calculation, the stray field of the Py film 
has been averaged over given thickness of NbSe2 film, 
and as such put into Eqs. (Q][2|). The simulation region 
was 2W x W large, and simulations were performed for 
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FIG. 4: (Color online) The equilibrium vortex phase diagram 
of the NbSe2 film at T — 4.2 K, as a function of thickness 
of the underlying Py film D and external perpendicular field 
H (the domain structure of the Py film is assumed to be 
unperturbed by the applied field). N denotes the number 
of vortex chains along the positive magnetic domains (illus- 
trated by |V>| 2 contourplots on the right, dark/light color - 
low/high density). At high fields, vortices penetrate areas 
above negative magnetic domains. 

f (0) = 10 nm and T = 4.2 0.6T C . The vortex states 
obey the imposed linear confinement and form multiple 
vortex chains (N of them) along the stripe domain with 
out-of-plane magnetization parallel to the applied field 
(as expected from the theory of magnetic pinning of vor- 
tices, see e.g. Ref. [24]). However, it should be empha- 
sized that in reality transitions between different TV-chain 
states are not of first order. In increasing field, added 
vortex causes local chain instability, by changing locally 
the transverse "optical" mode [25[. More instabilities 
eventually lead to a new ordered TV-chain state via a so- 
called "zig-zag" type transition characteristic for quasi- 
1D Wigner crystals [26j . These transverse instabilities 
leading to different ordering are evident in Fig. [2J where 
calculated vortex structures for the case of D = 1 /im 
are given for N = 1 — > 2 and N = 2 — > 3 chain rear- 
rangement, in comparison with experimentally obtained 
snapshots of vortex configurations. 

Note however that vortex spacing Ax along the chain 
prior to the N = 2 — > 3 transition remains larger than 
W/2(N + 1), a value predicted for N — > N + 1 transition 
in Ref. [9[. In Fig. ?? we show the calculated vortex 
spacing along the chain prior to the formation of a new 
chain as a function of the thickness of the magnetic film 
(i.e. proportional to the domain width). The main con- 
clusions following from Fig. ?? are that (i) threshold 
vortex spacing decreases with thickness of the magnetic 
film, but saturates for D > 1 /im, and (ii) Ax is very 
weakly dependent of the number of chains iV, contrary 
to findings of Ref. [9]. Note that in quasi- ID confined 
classical systems, with hard-wall lateral boundaries, Ax 
also decreases with N (see e.g. [27]), while in our system 
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lateral confinement is of parabolic type (where decrease 
of Ax with TV is slower [26j). 

The discrepancy from the model [9] is caused by 
the non-rigid lateral confinement of the vortex chains. 
Namely, the superposition of applied magnetic field and 
the stray field of the magnetic domains leads to slow ex- 
pansion of the confinement potential for the supercon- 
ducting vortices. From the magnetic force microscopy 
images at room temperature [8] and the STM maps in 
Fig. [2] we can reconstruct the position of the stripe do- 
mains under the NbSe2 layer, and the domain bound- 
aries are outlined with dashed lines in Figs. [2JSJ How- 
ever, beyond N = 3 one notices that vortex chains may 
already occupy areas slightly beyond the domain bound- 
aries. Therefore, with increasing magnetic field, not only 
that new vortices are added to the system, but confine- 
ment potential is broadened as well. This leads to further 
lateral instability of the chains, in the areas of strongest 
interaction with effective confinement. Moreover, the at- 
tractive force acting on newly added vortices may grad- 
ually become comparable to the overall long-range inter- 
action with earlier formed vortex chains, forcing the new 
vortices to reside in the anti-pinning channel above the 
magnetic domain of opposite polarity. This is shown in 
Fig. El for fields ±300 Oe in experiment and 325 Oe in the 
simulation. Due to competing interactions, several vor- 
tices are 'levitating' above the energetically non-favorable 
domain. Beyond the field of 325 Oe in the simulations, we 
were unable to recover the linear chain state, unless we 
had the commensurate number of vortices in each chain. 
In every other case, the attempt of the large number of 
vortices to form an Abrikosov lattice in competition with 
linear confinement, resulted in a sinusoidal-like modula- 
tion of the vortex structure, such as one shown in Fig. 

ETb). 

Further increase of the applied field leads to a larger 
number of vortices in the chains as well as above the 
anti-pinning domain. At certain conditions vortices at 
the anti-pinning channel may form a ID structure, be- 
ing equally repelled (on average) by chain structures 



from both adjacent domains. This is a very peculiar 
state, since vortex repulsion in the anti-pinning domain 
is strongest in its central region. The existence of such 
vortex ordering is verified in both theory (at w 385 Oe) 
and experiment (at 400 Oe), as shown in Fig.O Beyond 
this field, vortices have little difficulty trespassing to the 
anti-pinning domain, which finally leads to a somewhat 
distorted Abrikosov lattice across the whole film. 

To summarize, we studied the Abrikosov vortex 
configurations in magnetically coupled superconductor- 
ferromagnet bilayers. We elucidate the mechanism of 
quasi- ID topological vortex phase transitions in this sys- 
tem by theoretical simulations and STM vortex imaging. 
We demonstrate that transitions between ID ordered vor- 
tex phases are accompanied by local transverse instabil- 
ities resulting in local vortex disorder. This, in combi- 
nation with magnetic field-relaxed confinement, leads to 
lower threshold vortex densities for the geometrical tran- 
sitions than in strongly confined systems such as super- 
conducting stripes and hard- wall classical clusters. At 
higher magnetic fields the change of number of vortices 
is accompanied with the change in confinement potential 
resulting in nonlinear system response and vortex chain 
formation above the domains with opposite polarity. As 
a result, we found the 'levitation' of vortices above the 
magnetic domains of opposite polarity, as well as vortex 
chain formation in this energetically unfavorable region, 
while the earlier formed chains are being sinusoidally dis- 
torted from the equilibrium Abrikosov vortex lattice. 
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FIG. 5: (Color online) Average vortex spacing along the 
chain, prior to N — > N + 1 chains configurational transition 
(color coding corresponding to Fig. |4}, as a function of the 
thickness of the Py film. 
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